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Abstract 

Ten microsatellite loci were isolated and characterized for Stys's bush-cricket, Isophya stysi 
Cejchan (Orthoptera: Tettigoniidae), an endemic Orthoptera species to the Carpathian Basin, us- 
ing an enriched genomic library procedure. The polymorphism of these loci were tested in two 
populations of /. stysi, and the number of alleles per locus varied from 4 to 16. The expected and 
observed heterozygosities ranged from 0.612 to 0.925 and from 0.625 to 1.000, respectively. The 
interspecific applicability of these microsatellites was evaluated by amplification in 20 related 
species: Isophya camptoxypha, Isophya sicula, Isophya ciucasi, Isophya pienensis, Isophya harzi, 
Isophya kraussii, Isophya zubovskii, Isophya rectipennis, Isophya modesta, Isophya longicauda- 
ta, Isophya dobrogensis, Isophya hospodar, Isophya speciosa, Isophya modestior, Poecilimon 
fussii, Poecilimon affinis, Polysarcus denticauda, Barbitistes constrictus, Leptophyes discoidalis, 
Phaneroptera falcata. All primer pairs for the 10 loci yielded successful amplifications in at least 
one other taxon from the Isophya genus. This set of microsatellite loci would be useful for genet- 
ic studies in /. stysi and other species of the genus Isophya. 
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Introduction 



The genus Isophya (Orthoptera: Tettigoniidae) 
is one of the largest of the Palearctic Orthop- 
tera, comprising a large number of species 
with a high degree of morphological similarity 
(Heller 1988; Warchalowska-Sliwa et al. 
2008; Chobanov 2009; Grzywacz-Gibala et al. 
2010). The genus has a high number of en- 
demic species (Sevgili 2003; Sevgili and 
Heller 2003), most of them having restricted 
distribution ranges, limited mobility, and spe- 
cific topographic requirements (Sevgili et al. 
2006). Many European Isophya species are 
found in isolated populations with low densi- 
ties of individuals, displaying habitat 
dependence due to their feeding preferences 
for dicotyledonous plants (Bauer and 
Kenyeres 2006). These characteristics make 
them vulnerable to anthropogenic disturb- 
ances that lead to the reduction of their natural 
habitat, and many of the species have been 
rendered endangered. 

Isophya sty si Cejchan lives in small popula- 
tions throughout its range and has a 
fragmented distribution (Orci et al. 2005). It 
lives in mesophylic grasslands, near forests 
and forest clearings (Kis 1960; Pecsenye et al. 
2003), and adults can be found mainly on high 
herbaceous plants and small shrubs (Iorgu and 
Iorgu 2008). The species is endemic to the 
Carpathian Basin (Kis 1960; Nagy 2005), be- 
ing described from the Carpathian Mountains 
in Slovakia, and is also found in eastern Hun- 
gary, Poland, Romania, and Ukraine (Heller et 
al. 2004). In Romania, it is common in the 
Apuseni Mountains and Transylvania, and 
rare in the Oriental Carpathians and in the for- 
ests of the Moldavian Plateau (Iorgu et al. 
2008). 
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/. stysi is protected by national and European 
laws (present on Annex II of EU Habitat Di- 
rective) and requires strict conservation 
measures. In order to ensure effective conser- 
vation management, the genetic diversity of 
the populations of this endangered species 
needs to be evaluated. Microsatellite DNA is 
an optimal molecular marker for studies of 
genetic diversity in natural populations as it 
can assess population structure (Goldstein and 
Schlotterer 1999). Population genetic analyses 
can also be used to identify management units 
based on ecological and genetic variation, and 
to trace threatened populations in need of con- 
servation priority (Wan et al. 2004; Palsboll et 
al. 2007). The major aim of this study was to 
describe the first microsatellite loci for /. stysi 
and to report the results of cross-species am- 
plification tests in 20 other related Orthoptera 
species. 

Materials and Methods 

The isolation of the microsatellite loci for /. 
stysi was performed following a standard pro- 
tocol for the construction of a microsatellite- 
enriched library (Bloor et al. 2001). Genomic 
DNA was isolated from the hind femurs of 
two individuals of /. stysi using a phenol- 
clorophorm protocol (Sokolov 2000). 

Approximately 10 ug of genomic DNA was 
digested using Sau3A\ restriction enzyme 
(Fermentas UAB, www . fermentas . It ) . Adap- 
tor-ligated DNA fragments ranging from 400 
to 1000 bp were selected, and enriched using 
3' biotin-labelled CA and GA repeat oligos 
bound to streptavidin coated magnetic beads 
(M-280 Dyneabeads, Dynal, Invitrogen, 
www.invitrogen.com) . The DNA fragments 
were then ligated into the p JET 1.2 vector 
(Fermentas UAB) and transformed into DH5a 
Escherichia coli competent cells for cloning. 
The enriched genomic library was screened 
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for repetitive sequences, and 73 clones con- 
taining inserts with microsatellite motives 
were selected. These were further sequenced 
using the LICOR 4300L Genetic Analyzer. 
The similarity of the flanking regions and the 
microsatellite length were determined with 
SciRoKo 3.4 (Kofier et al. 2007), and 25 se- 
quences were chosen suitable for primer 
design. Primers were designed using Primer 3 
program (Rozen and Skaletsky 2000). 

All 25 primer pairs were initially tested for 
amplification using 8 individuals of /. stysi. 
Primers that yielded products of expected siz- 
es were given an Ml 3 sequence tail to allow 
analysis on the LICOR 4300L genetic analyz- 
er (M13F: 5'-cacgacgttgtaaaacgac-3', M13R: 
5 '-ggataacaatttcacacagg-3 '). 

The genotyping reactions were performed in a 
10 uL reaction volume, containing about 30 
ng of DNA template, 10 mM Tris-HCl (pH 
8.8 at 25° C), 50 mM KC1, 0.08% (v/v) Non- 
idet P40, 1.5 to 2.5 mM MgCl 2 (see Table 1 
for details for each locus), each dNTP at 0.1 
mM, each primer at 0.1 uM, 0.02 uM of 
IRD700, or IRD800 labeled Ml 3 primers (the 
same sequence as the M13 tails), and 0.5 units 
of Taq DNA polymerase (Fermentas UAB). 

The PCR program used consisted of 5 minutes 
denaturating at 95° C, followed by 30 cycles 
of 30 sec at 95° C, 30 sec at the annealing 
temperature (see Table 1 for each locus), and 
45 sec at 72° C, ending with a 7 min final 
elongation stage at 72° C. The genotyping 
process was performed using the Saga GT 3.2 
software package (LI-COR Biosciences, 
www.licor.com) . 

The degree of polymorphism of the 10 select- 
ed loci was tested in two populations of /. 
stysi collected from two mesophylic meadows 
in Romania (23 individuals from Ceahlau 
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Mountain, Neamt, County: 47° 01' 14' 'N, 25° 
57' 16" E; 21 individuals from Nucsoara, Alba 
County: 45° 29' 09" N, 22° 56' 03" E). Ge- 
nomic DNA was extracted from the middle 
leg of each individual using the Nucleospin 
Tissue kit (Macherey-Nagel, www.mn- 
net.com) . 

The null alleles frequencies were estimated by 
a maximum likelihood algorithm as imple- 
mented in FreeNA (Chapuis and Estoup 2004; 
Chapuis et al 2008), and tests for linkage dis- 
equilibrium were carried out using GenePop 
v. 4.0.10 (Raymond and Rousset 1995; Rous- 
set 2008). 

Results and Discussion 

Only ten out of 25 primer pairs proved to be 
polymorphic and were deemed acceptable for 
population genetic studies. The number of al- 
leles at each polymorphic locus, their size 
ranges, observed and expected heterozygosi- 
ties, as well as deviation from the Hardy- 
Weinberg equilibrium were calculated using 
GenAlEx 6.4 (Peakall and Smouse 2006) and 
are summarized in Table 1 . The microsatellite 
loci showed high levels of polymorphism, the 
number of alleles per locus ranging from 4 to 
13 in the population from Ceahlau Mountain 
and from 4 to 16 in the population from Nuc- 
soa-ra. In the population from Ceahlau Moun- 
tain, the observed and expected 
heterozygosities ranged from 0.625 to 0.957 
and from 0.612 to 0.891 respectively, with an 
average of 0.843 and 0.799 respectively. In 
the population from Nuc-soa-ra the observed 
and expected heterozygosities ranged from 
0.714 to 1 and from 0.684 to 0.925 respective- 
ly, with an average of 0.876 and 0.844 
respectively. 

Significant deviation from the Hardy- 
Weinberg equilibrium was observed in 3 out 
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of 20 possible single exact locus tests (p < 
0.05), IST3 only in the Ceahlau Mountain 
population, and IST18 and 1ST 23 only in the 
Nuc-soa-ra population. Null alleles were esti- 
mated as present in IST15 (estimated 
frequency f.e. = 0.028), IST21 (f.e. = 0.059), 
and IST24 (f.e. = 0.019) loci in the Ceahlau 
Mountain population, and in IST23 locus (f.e. 
= 0.073) in the Nuc-soa-ra population. These 
results, together with the relative small sample 
size (Ceahlaua: 23; Nuc-soa-ra: 21), may ex- 
plain the deviation observed in some of the 
Hardy- Weinberg equilibrium tests. No signifi- 
cant linkage disequilibrium was found 
between loci pairs in tests performed across 
all populations. 

The molecular variance analysis, calculated 
using GenAlEx 6.4, showed significant differ- 
entiation between the two populations (p = 
0.01), with a moderate pairwise Fst value of 

0. 056. The genetic differentiation between the 
two populations can be explained by the geo- 
graphical distance between them (almost 290 
km in a straight line), which can determine a 
low level of geneflow. 

In order to assess interspecific amplification, 
the polymorphic loci were also tested in 20 
additional species from the Phaeropteridae 
family: 14 species of the genus Isophya (I. 
camptoxypha, I. sicula, I. ciucasi, I. pienensis, 

1. harzi, I. kraussii, I. zubovskii, I. rectipennis, 
I. modesta, I. longicaudata, I. dobrogensis, I. 
hospodar, I. speciosa, I. modestior), two spe- 
cies of the genus Poecilimon (P. fussii, P. 
affinis), and four species from different genera 
(Barbitistes constrictus, Polysarcus denticau- 
dus, Leptophyes discoidalis, and 
Phaneroptera nana) (Table 2). DNA samples 
from two individuals of each species were 
genotyped using the same PCR conditions 
used for /. stysi. All primer pairs amplified in 
at least one other taxon from the Isophya ge- 
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nus, and only 6 of them amplified for species 
outside of the genus (IST2, IST5, IST6, IST9, 
1ST 15, and IST21). IST2 and IST5 loci ampli- 
fied in all Isophya species tested and 1ST 15 
locus amplified in 12 species. Eight of the mi- 
crosatellite loci (Table 2) amplified in /. 
modestior, which is considered closely related 
to /. stysi from morphological and bioacustical 
data (Warchalowska-Sliwa et al. 2008). 

These data show that the microsatellite mark- 
ers isolated for /. stysi may prove to be very 
useful in population genetic studies on other 
species of the genus Isophya, but their poten- 
tial for cross-species amplification is limited 
outside the genus. These novel polymorphic 
loci should be a useful tool to study the genet- 
ic diversity and structure of I. stysi 
populations and to develop better conservation 
measures for this endangered species. 
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Table I. Characterization often microsatellite loci developed for Isophya stysi. 



Locus 


Gen Bank Acc. No. 


Primer sequence 


Repeat motir 


In (°C) 


Mgl"l; (lIlM) 


N 


Na 


Size- range 


Ho 


He 


HV\ F.-p 


IST2 


JQ41323I 


F • AC TT AT( " T( ,CC T( *C(\C,C,TCT 

1 . / \ V 1 1 / \ 1 V- 1 V1V.. V- 1 V- V. V IV IV J 1 V 1 


(AC),, 


54 


1.5 


36 


6 


176-188 


0.78 


0.67 


0.464 


li- CCCTCA ACTCCA ACC, A ATfiT* 
i\. v J v v_ i v. rtrtvj 1 v J v r\r\vjvjr\r\ i vj I 


1ST3 


JQ4 13235 


T: TATGGCAAACCCCGATATTT* 


(AG)„TG(AG),,G 
G(AG), 


51 


1.5 


36 


18 


145-251 


0.81 


0.91 


0.008 


R: TTGTAGCC 1 CAGTGGCATGT 


IST5 


JQ4 13236 


1 : CAACGCATGGGACACACTAT 


(CT),CGGT(CT)., 
CC'(CT) 4 


53 


2.5 


36 


13 


228-282 


0.83 


0.83 


0.102 


R: TGCCCGTCCATAACACACTA* 


1ST6 


JQ4 13237 


F: GGTACTCCTGCAGGCTGAAA 


(CT).„ 


51 


2.5 


36 


15 


148-228 


0.89 


0.9 


0.159 


R: CAGAAGCATAATCCATCACAGG' 


IST9 


J04 13238 


F: TTTAACGGGCATCAGGTTTC 


(AC)„GC(AC), 3 


53 


2 


35 


14 


219-253 


0.86 


0.9 


0.377 


R: CAAAGTCGACGCAGAACAAA* 


IST15 


JQ4 13229 


F: ACGGATGGATTGAATGAAGC 


(CT),CACT 

1 V , \ J i i M v_ i \ f 


51 


2 


35 


18 


204-284 


0.91 


0.91 


0.191 


R: CCAACCAAACGACACCACTA 


ISTI8 


JQ4 13230 


F: ACGAACAAACGAACGAAAGC 


(TG)„CG(TG), 0 


50 


2.5 


36 


10 


208-286 


0.86 


0.78 


0.012 


R: AAGGCCGCATITCTACTTCA 


IST21 


JQ4 13232 


1 : CCAAACAGAAGACAGCTCG IT 


(CA) 17 


55 


2 


33 


IS 


235-289 


0.82 


0.88 


0.147 


R: ACCCACCACGAGGTAGACAG* 


IST23 


JQ4 13233 


F: GATTCAGGCAAGCGGTAGTC" 


(TC),„CC(TC)„ 


54 


1.5 


36 


15 


145-205 


0.81 


0.89 


0 


R: GCTTiflTGTGAf'GAAACCAC 


IST24 


JQ4 13234 


F: ACTTTAAGTTTCGCCGAGCA 


(AC)v.(AG)u 


53 


2 


36 


21 


165-261 


0.83 


0.92 


0.026 


R: ATGAGTTCCGTCCGTCAGTC* 


Mean 














14.5 




0.84 


0.86 




SF 














1.34 




0.01 


0.03 





*: specifies M 1 3R tailed primer; # : specifies M 1 3F tailed primer; Ta: annealing temperature (° C); MgCb: MgCb concentration in the 
PCR reaction (mM); N: number of tested individuals; Na: number of alleles; Ho: observed heterozygosity; He: expected heterozy- 
gosity; HWE-p: uncorrected exact f>-value of the Hardy-Weinberg equilibrium test. 

^ ' 



Table 2. Cross-amplification of the microsatellite loci in 20 species of the family Phaneropteridae. 



TAXON 


1ST 2 


1ST J 


1ST 5 


1ST 6 


1ST ') 


1ST 15 


1ST 18 


1ST 21 


1ST 23 


1ST 24 


Vi 


sr 


Na 


sr 


Na 


sr 


Na 


sr 


Na 


sr 


Na 


sr 


Na 


sr 


Na 


sr 


Na 


sr 


Na 


sr 


Isophya contptoxypha 


4 


166-182 




147-175 


4 


206-236 


4 


158-218 




215-247 






1 


217 


2 


253-293 


2 


137-141 






Isophya sicula 


4 


170-180 


2 


157-185 


2 


202-238 


3 


184-236 






1 


240 










2 


139-153 






Isophya ctliCCtsl 


4 


138-182 


2 


157-169 


3 


210-236 


4 


156-226 


2 


255-247 


3 


146-272 


1 


219 


4 


249-279 










hophyaptenensts 


3 


176-182 


4 


145-189 


3 


198-282 


3 


204-220 






3 


224-242 


1 


217 






2 


141-147 






Isophya harzi 


2 


168-180 








206-236 






4 


223-247 


2 


232-250 


3 


203-221 


4 


1 29-295 


2 


151-159 






Isophya kraussi 


2 


170-180 






2 


232-236 


4 


180-234 






4 


214-242 










4 


133-163 






Isophya -iibovskii 


2 


178-182 


3 


175-191 


2 


202-236 


4 


2 1 6-234 






4 


232-256 




219-297 














Isophya rectipennis 


2 


174-182 






2 


212-236 


4 


174-230 


1 


147 


2 


220-230 


2 


221-295 














Isophya mojesta 


4 


178-194 


2 


153-159 


1 


284 






2 


213-247 


3 


236-262 


I 


221 


2 


221-267 


2 


147-59 






Isophya longicaudata 


2 


176-182 


3 


155-177 


! 


236 


4 


158-238 






4 


252-272 


1 


221 


4 


223-271 


1 


163 






Isophya dobrogensis 


2 


176-182 


3 


127-179 


2 


260-274 






3 


205-219 


2 


242-260 


1 


223 


3 


237-295 










Isophya hospodar 


3 


168-182 






1 


294 


1 


212 










2 


163-221 


4 


253-287 










Isophya speciosa 




164-180 


2 


157-175 


2 


150-166 


3 


242-252 






3 


226-238 


2 


131-223 


2 


137-293 


4 


139-183 






Isophya modestior 


2 


176-182 


4 


153-197 


3 


262-288 






4 


199-227 


2 


212-240 






4 


133-285 


4 


149-177 


4 


189-263 


Poecilimon fussi 




176-186 








152 






1 


247 






















Poeciliinon affinis 


1 


158 








254 


1 


298 


























Barhit isles constriclits 






























1 


253 










Polysarciis denticaudus 


















1 


219 










2 


183-255 










Leptophyes discoidaiis 












250 


2 


210-214 


























Phaneroptera nana 












394 












240-278 



















Na: number of alleles; sr: size-range; -: no amplification or multiple bands. 
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